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Abstract 



A method of MC simulations including quantum interference, pro- 
\Q ' posed recently by A.Krzywicki and the present author, is explained. 

On' 

1. Introduction. 

The aim of this talk is to explain details of a scheme for MC simulations 
of multiparticle production including HBT interference which we proposed 
recently together with Andrzej Krzywicki [|lj] . The problem became suddenly 
of great practical importance, when it was realized that the HBT effects may 
seriously affect some precise measurements of the standard model parame- 
ters @. As we have heard yesterday from Krzysztof Fialkowski, the existing 
implementations of quantum interference into standard MC codes suffer from 
many problems, theoretical as well as practical ones |3], |J. Therefore con- 
struction of a viable MC code correcly including the interference effects is 
badly needed. To do this, however, it seems first necessary to formulate the 
problem and the goals to be achieved. It should perhaps be emphasized at 
this point that the effects of HBT interference are by no means "trivial" 
or "automatic", as it is sometimes believed. On the contrary, they depend 
in essential way on the physics of the problem. Consequently, there is no 
single, unique method of implementing the HBT interference into existing 
codes. One must therefore be careful to spell out the underlying physical 
assumptions. 
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We begin in the next section by a brief reminder of the relation between 
the distributions in x-space and in momentum space which will allow to intro- 
duce the necessary concepts and to formulate the assumptions. In Section 
3 the physical meaning of the procedure is explained in terms of Wigner 
functions. Some comments and outlook are given in the last section. 

2. Density matrix and relation between x-space and p-space. 

Let 4>(qi,q2, ■■■Qn,C() = i>(q,(x) be the probability amplitude for produc- 
tion of N particles with momenta qi, q%, ....qN = q. a denotes a collection of 
all other quantum numbers which may be relevant to the process in question 
(they may be, e.g., the momenta of other particles which we do not wish to 
consider explicitly in a "semi-inclusive" measurement). The density matrix 
in momentum space is then 



This matrix gives all available information about the system in question. The 
observed spectrum of particles reads 



We see from this formula that measurement of the momentum spectrum pro- 
vides only a rather limited information about the system: only the diagonal 
elements of the density matrix are determined. 

Let us now consider the coordinate space. We write 



where ip(x,a) is the probability amplitude for producing N particles at the 
points (xi, X2, ...Xn = x) and < q | x > is the known transformation ma- 
trix between momentum and coordinate space which shall be specified later. 
Introducing the density matrix in coordinate space 




(1) 






(3) 





we obtain the relation 



p(q, q') = J 



dxdx' < q | x > p(x, x ) < x \ q > 



(5) 
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which shows that transformation between description of the system in mo- 
mentum and in coordinate space requires the knowledge of the full density 
matrix. The measured distributions (which give only the diagonal elements) 
are not enough. 

To continue, we need an explicit form of the transformation matrix < q | 
x >. As it is different for identical and non-identical particles, we shall treat 
these two cases separately. 

(i) non-identical particles 

In this case 

< q | x >= expiiqx) = exp[i(qiXi + g2^2 + ••■?atXjv)] (6) 

(all powers of 2ir are included in normalization of dx and dq). Substituting 
this into @ we have 

p (q,q') = J dxdx'e t{qx - q ' x,) p(x,x') = J dx + dx~ e l{q ~ x+ +q+ x ~ ] p{x , x') (7) 

where 

q + = ~(q + q');q~ — q — q'] x + = ~(x + x') ; x~ = x - x' (8) 

(from now on we denote the quantities referring to non-identical particles by 
a subscript 0). 

From (0) we obtain for the spectrum of non-identical particles 

^o(?) = Po(q,q) = J dx + dx~e iqx ~ p(x + ,x~). (9) 

This formula shows explicitly that the measured momentum spectrum of 
non-identical particles does not give any information on distribution of par- 
ticles in coordinate space: the x + dependence is integrated over. Instead, 
we obtain information on x~ dependence, i.e., to what degree the density 
matrix in coordinate space is non-diagonal. In more physical terms, the mo- 
mentum spectrum gives information only on the coherence properties of the 
system in the coordinate space. Indeed, the off-diagonal part of the density 
matrix measures how much the result of the integration over the internal 
quantum numbers of the system (denoted by a in the Eq.(^)) is affected by 
the cancellations due to "incoherent" summation of terms with "randomly" 
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distributed phases. For illustration, let us consider a simple parametrizaton 
of the density matrix in the form 

p(x + , x~) = p{x + )e^ x ~ )2/l ' {x+) (10) 

where p(x + ) is the distribution of particles in coordinate space and l c (x + ) is 
the "coherence lenght". In this case we obtain 

n (q) = J dx + p(x + )ll(x + )e- Aq2l ^ x+) (11) 

and we see explicitly that the momentum spectrum measures the average 
value of a quantity depending on the coherence lenght l c (x + ). 

This point is dramatically emphasized for a limiting case of the system 
which is fully incoherent in coordinate space, i.e. for which the density matrix 
is purely diagonal (l c — » 0). In this case we obtain the spectrum which is 
entirely independent of q! (see e.g. || for a more detailed discussion of this 
result and other effects of incoherence). 

One final remark about normalization: one sees from (H) that 

/ n ( q)dq = J dx+p(* + , x- = 0) = / dxp(x, x). (12) 
(ii) identical particles. 

In this case we have to symmetrize the transformation matrix over particle 
momenta and positions and thus we obtain 

(MpX> i?PS ( 13 ) 

where P is a permutation of the numbers (1,2,. ..N) and qp are the momenta 
(qi, q2, (/at) ordered according to the permutation P. Introducing ([13]) into 
(fD we have 




Using (0) this can be rewritten as 

p(q, q') = ^Y1 Po(gp, q'p>) (15) 
iV • p,p' 
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so that we obtain for the momentum distribution of the identical particles 

Q(q) = p( g> q ) = 2_ £ /" dx + dx~e i{q pp> x++q pp> x ^p( x + , x~) (16) 
^ * ■ p pi j 

where qp~ P , = \{qp + qp>) and qp P , = qp — qp>. Eg . (|T6|) shows explicitly that 
momentum distribution of identical particles gives information on both x~ 
and x + dependence of the density matrix. For the example (pTOj ) one obtains 

Sl(q) = E / ^ + e^' x+ p(^ + )^ 3 (^ + )e~ 4{ ^' )2 ' c2(x+) (17) 

which clearly shows that the dependence on momentum differences qp P , is 
sensitive to x + dependence of the particle density in coordinate space p(x + ) 
and of the "coherence length" l c (x + ). 
Three remarks are in order. 

(a) The momentum spectrum of N particles given by ( ]T5| ) is expressed in 
terms of the density matrix of N particles in the coordinate space and thus 
cannot be reduced (without further assumptions) to the expression involv- 
ing only single particle density. In particular, it depends on all N-particle 
correlations in the coordinate space. Usually these correlations are neglected 
(i.e. the density matrix p(x + ,x~) is written as a product of single particle 
matrices). Although this is a reasonable procedure in the absence of any adi- 
tional information, it should be kept in mind that future data may reguire 
to include these correlations Jl|, 0] . 

(b) The normalization of the spectrum ([16]) is different from that of non- 
identical particles given by ([12]). Integration over particle momenta gives 




(18) 

This result shows that the guantum interference changes not only the dis- 
tribution of produced particles but also the production cross-section (i.e. it 
acts as final-state interaction). 

(c) When all particle momenta are egual to each other we obtain from 
© and (0) 

Q(q) = N\n (q) if qi = q 2 = ... = q N (19) 
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consistent with the standard treatment. 
3.Wigner functions. 

The density matrix has a clear physical meaning, as seen from Eqs. (j^) 
and (|4]). Its intuitive meaning is, however, more difficult to grasp. Therefore 
it is useful to consider a Fourier transform 

W(q + ,x + ) = J dx~e iq+x ~ p(x + ,x-) (20) 

which is the generalization of the well known Wigner function (defined usu- 
ally for single particle spectrum). It is seen from ( p0|) that W(q + , x + ) is a 
quantum-mechanical generalization of the classical particle density in mo- 
mentum and in coordinate space (Boltzmann phase-space density). 

Using ( p0|) and ©,(0) the particle densities for non-identical and iden- 
tical particles can be respectively written as 

n Q (q) = J dxW(q,x) (21) 

n (ff) = M E / dxW(^±^,x)e^-^ (22) 

From these relations one sees again explicitly that while the momentum dis- 
tribution of non-identical particles does not give any information on the par- 
ticle distribution in x-space, the measured momentum spectrum of identi- 
cal particles is sensitive to x-dependence of the Wigner function, i.e. to 
x-dependence of the distribution. 

The advantage of using the Wigner functions is that they appeal to one's 
intuition (being the analog of the Boltzmann distribution) and thus the re- 
sulting formulas are easier to interpret. Of course it should be kept in mind 
that this analogy is limited by the fact that a Wigner function is -in general- 



locally not positive definite. It can oscillate and, as seen from (21), the oscil- 
lations cancel out only after integration over x. However, these oscillations 
can play a significant role in the Eq.(|22|) for identical particles by conspir- 
ing with oscillating terms in the integrand to contribute significantly to the 
result. This is how quantum mechanics shows up in the problem. Thus re- 
garding Wigner functions as Boltzmann phase-space density is possible only 
when they are appropriately smoothed out to remove the oscillations. The 
price to pay is that, in general, the resulting probabilistic description can 
only be trusted for when the momentum differences in ([22]) are not too large. 
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4. The proposal for MC simulation. 

Standard Monte Carlo algorithms generate multiparticle events according 
to an assumed model for the momentum spectrum Qo(q) which does not 
include quantum interference. The problem is to correct the weights of these 
Monte Carlo events once they were generated. 

It is clear that this cannot be achieved without additional assumptions. 
Our proposal is to assume that the corrected spectrum is given by Q(q) of 
the Eq.(|22|) with the same Wigner function as that present in fl2T|) . 

I would like to emphasize that this assumption is far from obvious, al- 
though it is usually accepted without further comments (see, e.g.,||). It 
assumes that the identical and non-identical particles are produced in the 
same way. Clearly, this can only be an approximation (resonance production, 
for example, influences differently identical and non-didentical particles). As 
discussed by Bo Andersson at this meeting, it is also violated -generally- in 
the Lund model |J. Hopefully it is not unreasonable for events with many 
particle which we are concerned with Q. 

For an effective use of the Eqs. ( |2l|j2^ ) we need an expression for the 
Wigner function which reproduces the spectrum fio(<z) f° r non-identical par- 
ticles. Therefore we write 

W(q,x)=n (q)w(q,x). (23) 
It follows from ([H]) that w(q, x) obeys the normalization condition 

w{q, x)dx = 1. (24) 



We see that w(q,x) is the quantum analog of the conditional probability: 
given that particles with momenta q\, q%, qN are present in the final state, 
w is the probability that they were emitted at the points X\,X2, ...,xn- 
When (^) is inserted into (|22|) we obtain for the correcting weights 

s{q) =m = ^^^r wi ^' qp ~ qp,) (25) 

where 

w(q, A) = J dxw(q,x)e iAx (26) 

lr To avoid this assumption one needs either a specific model of multiparicle amplitudes 
(see e.g. |6|) or a direct calculation from the first principles. 
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with 

w(q,0) = l. (27) 

Clearly, w(q, x) is rather arbitrary and must eventually be determined by 
analysis of the data. In absence of any information, and to exploit fully the 
intuitive meaning of Wigner functions, we propose -as a first step- to neglect 
possible oscillations and to take w(q, x) in a form which is everywhere positive 
definite, so that it can indeed interpreted as a probability distribution. 

The formula (|25|) cannot be used at it stands for most of the existing MC 
algorithms because they use an iterative procedure which provides Qo(q) 
only for one given set of momenta and not all the sets needed in (|25]). This 
difficulty can be dealt with by observing that one does not make a big error 
by replacing in (|25|) £l ( qp+ 2 qp ' ) by Q (qp). Indeed, those terms in (|25|) where 



this approximation is poor are suppressed by the rapidly decreasing factors w 
and thus need not be calculated with great precision. Eq.(p5|) now becomes 



^■sj*^'"-""' (28) 

The same argument can be used to see that the weights given by ([28|) 
are positive, as required for MC simulations. To this end we observe that, as 
seen from ([27]), they are certainly positive if the difference between particle 
momenta are small. Thus the positivity is quaranteed in the region where 
our approximation for Wigner functions is valid. As we have argued before, 
outside of this region the non-diagonal w( qp £ p ' , qp — qp>) are small and 
-whether positive or not- do not play any role in the sum (f28|). 



To proceed, further working assumptions are needed. In |l] we proposed 
to start with w (q, x) in a factorized form, each factor being a superposition 
of exponentials. Such factorization is most likely not exact (some indications 
of this were shown by Hans Eggers at this meeting |7j) and may be corrected 
when the actual data are fitted. For more details we refer the reader to 
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